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A METHOD FOR THE PROCESSING AND ANALYSIS 


OF DIGITAL TERRAIN ELEVATION DATA 
By Bobby G. Junkin* 

SUMMARY 

This report presents a method for the processing and 

analysis of ^ digital topography data that can subsequently 

» \ 

be entered in an interactive data base in the fora of slope, 
slope length, elevation and aspect angle. Included are a 
discussion of the data source and specific descriptions of 
the data processing software programs. In addition, the 
mathematical considerations involved in the registration of 
raw digitized coordinate points to the UTM coordinate system 
are presented. Scale factor considerations are also included. 
Results of the processing and analysis are illustrated using 
the Shiprock and Gallup Quadrangle test data. 


*Earth Resources Laboratory, NASA/NSTL, Slidell, Louisiana 70458 
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INTRODUCTION 


Land resource managers are becoming increasingly aware 
of the problem of converting disparate sources of data in map 
format into a form suitable for processing on a computer- 
oriented information system. This information, acquired from 
map sources or remote sensor data obtained from aircraft and 
satellites, in compiled into data bases which contain infor- 
mation on land use, topography, soil, rainfall, population 
density, etc. This report defines the procedures and tech- 
niques in use at the NASA/ERL for processing digital topography 
data that can subsequently be entered in a data base in the 
form of slope, slope length, elevation and aspect. 


DIGITAL TERRAIN ELEVATION DATA 


The National Cartographic Information Center (NCIC) of 
the U. S. Geological Survey, Department of the Interior, provides 
a national information service to make cartographic data of the 
United States available to the public and to various federal, 
state, and local agencies (reference 1). These cartographic 
data include standard 9 -track 800 BPI or 1600 BPI computer- 


compatible tapes (CCT) which contain digital representations 
of terrain elevations. These tapes, which are produced by 
the Defense Mapping Agency Topographic Center (DMATC) from 
the 1 : 250 ,00C-scale series of maps, are copied and distributed 
to users by rhe NCIC. 

The DMATC utilizes a processing system which collects 
data from a 1 : 250 ,000-scale map using digital graphic recorders 
tied into a central processor with disk storage, magnetic 
tape output and verification plot capability. The function 
of this sytem is to generate a magnetic tape file containing 
a matrix of elevation readings extracted at 250 micrometer 
intervals. The 1 :250 ,000-scale map generally covers one 
by two degrees of longitude and latitude. The DMATC prepares 
two one-degree by one-degree matrices for each quadrangle. 

Each block of data is stored on a 7- track UNIVAC 1108 
computer tape by the DMATC and forwarded to the NCIC. The 
NCIC subsequently takes these data and stores up to eight 
one- degree quadrangles of longitude and latitude on each 
9 -track (1600 BP I) tape (four using 800 BPI tapes) . This 
latter data format is generally used by the ERL as its source 
of topographic information. 

General information concerning the source of data is 
given in Appendix D. 
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DATA PROCESSING SOFTWARE DEVELOPMENT 


General Program Functions 

The NASA/ERL software programs for geographical data 
analysis and display consist of four separate programs for 
processing digital terrain data tapes produced by the DMATC. 

These programs perform the following functions: (1) Transforms 

data in the local (x,y) digitizer system to the (Xg,Y N ) UTM 
system, (2) translates the 1° x 1° data set to an input origin, 

(3) adjusts the input data set to any desired output cell size, 

(4) computes slope, slope length, and aspect for each cell by 
use of the maximum gradient from a center cell to the surrounding 
eight cells, (5) generates output data files (elevation, slope, 
slope length, and aspect) which are in a form suitable for 
entering into the ERL Interactive Data Base Display Program (IDDP) . 

Specific Program Descriptions 

The processing of digital terrain tapes requires the 
execution of four separate programs. These programs are described 
herein. 

Program TQPREF . This program prints out pertinent infor- 
mation for each column of data as written to the reformatted 
output tape. The starting x-y points, ending y-point, number 
of y-points , and the first three and last three elevation data 
values are printed for each x-column. Also, the minimum and 
maximum elevation values are printed. This program reads the 
NCIC data tapes of 15840-word records (16-bit words) and writes 
a tape file of 2005-word records in the following format: 

t 
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word 1=0 

word 2 = record number 

word 3 = x coordinate in .01 inches for this record 
word 4 = starting y coordinate in ,01 inches 
word 5 = N = number of elevation values in this record 
words 6 through (6 + N) = elevation values of points 
(x,yi> through (x,y N > 

words (6 -4- N + 1) through 2005 = zero fill 

Program TOPTWQ . This program prints out all input control 
parameters which define the area of interest to process, column 
and row bias's, etc. Also printed out is a table of computed 
easting and northing and elevations for the first and last data 
point for each row. The reformatted tape from TOPREF is read. 
Card input control parameters are used to extract only that 
data required for processing the selected area and writes out a 
file suitable for sorting. Output from this program consists of 
7200-word records in 3-word groups of Xg and Y^ coordinate values 
for each cell and its elevation value . 

Program TQPSRT . This program prints out check point row 

numbers for raw data that are sorted and output in reverse order. 
The data that are written in T0PTW0 are sorted so that all 
elevation values corresponding to Yjj are written as the first 
record, those elevation values corresponding to Yjg + ^ as record 
2, etc., continuing on until record N = all elevation values 
corresponding to Y^ . These data are written to an intermediate 
scratch disk file, and then output to tape. 
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Program TOPODB . This program prints input control parameters 
and input levels for elevation and slope. Optionally, a printer 
map is printed showing coded levels for each data element in 
row- column format for all four types of output. Output consists 
of four files of data to disk for input to the IDDP: 

1. Elevation data 

2. Slope data 

3. Aspect data 

4. Slope length data 


MATHEMATICAL ANALYSIS 


UTM. Coordinate Equations 


Consider the example of a 1 : 250 ,000-scale map as shown 
in figure 1. The UTM (X £ ,Y N ) grid system on these type maps 
is not normally aligned with the (x,y) coordinate system of 
the digitizer system. Thus , an angle correction is required 
for each 1° x 1° area or file of data. The sheet corners indi- 
cated by the arrows are translated into the origin of the data 
file by the following: 


~ f 



x i = yi - y 9R i 


I = 1, 2 


( 1 ) 


where : 


(Xi , y^) = plate coordinates of digitizer 
points , in inches 
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ORIGINAL PAQS $ 

OP POOR QUALITY 


Figure 1 - Exaaple of USGS 1:250,000 Scale 





(Xqri , ygj^) * coordinates of SW sheet comer 
of file 1, in inches 

(x 0 R 2 , y 9R2 ) = coordinates of SE' sheet corner 
of file 2 , in inches 

The (x^, y|) data are corrected for the A angle between the 
polyconic coordinate system and the UTM coordinate system to 
yield corrected (xY, yV) digitizer plate points. 

The relationship between the corrected digitizer plate 

points and the corresponding UTM coordinates are given by the 

\ 

following transformation: 


X Ei _ 


s Xi 

+ X 0RJ 

(2) 

Y Ni " 

y" 

s Yi 

+ Y 6RJ 

(3) 

where J = 1 , 

2, 

and: 




( x 0RJ* ^SRJ^ = coordinates of SW or SE sheet comer 

Sy^ = X scale factor, file 1 or 2 

Sy^ = Y scale factor, 1 file 1 or 2 

It should be pointed out that the above equations, in effect, 
register the plate coordinates of the left side 1° x 1° area 
of a quad map r.o the SE sheet comer and the subsequent UTM 
coordinates are merely determined with respect to the (X 0R 2> Y 0R2^ 
origin of the SE sheet comer. 


Derivation of Equations for Data Registration 
The (x,y) digitizer plate coordinate points are registered 
to the UTM coordinate system through the utilization of the A 
angle between the poly conic coordinate system and the UTM 
coordinate system. Consider first the geometry in figures 2 
and 3. These figures depect the location situations that are 
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Polyconic 


IX 



Figure 3. West of Zone CM 



possible regarding a quad sheet relative to the central meridian 
of a zone. Figures 4, 5, 6, and 7 show the geometry between the 
(x,y) coordinate system and the UTM coordinate system. 

We first consider the geometry in figure 4 as relates to the 
right side of a quad located east of the zone CM. From this 
figure : 


y" - + hi (4) 

and: 

x" = a^ cosA (5) 

a^ = a;, sin A (6) 

from AT]_: 

y ' = b^ cosA ' (7) 

b 3 = sinA (8) 

solving (7) for b^: 

b-^ = y' / cosA (9) 

substituting (9) in (8) : 

bg = y* sinA (10 

cosA 


also: 

a 4 = x' - b 3 (11) 

substituting (10) in (11) : 

a 4 = x' - y ' sinA (12) 

cosA 
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and (12) in (5) : 


x" = x' cosA - y' sinA (13) 

substituting (12) in (6) and the result in (4) and also (9) in (4) 
yields : 

y" = x' sin - y ' sin^A + y* (14) 

cosA cosA 


where the (x,y) points are translated to the (x 0 ,y 0 ) origin 
point by: 


X ' = X - X, 


= y _ 


(15) 


Thus, equations (13) and (14) give the registration of the (x,y) 
plate coordinate point relative to its (x",y") position in the 
UTM coordinate system. These equations are valid for the right 
side of a quad map and east of the zone CM. 


Consideration of the geometry in figure 5 yields the following 
equations for the left side of a quad map, east of zone CM : 

x" = x ! cosA + y' sinA ] 

o (16) 

y" = -x sinA - y 1 sin^A - f y» l 

cosA cosA ] 

These equations are also valid for the right side of a quad, 
west of the zone CM. This can be verified by consideration of 
the geometry in figure 6. Consideration of the geometry in 
figure 7 also shows equations (13) and (14) to be valid for 
the left side of a quad map, west of the zone CM. 
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Scale Factor Equations 


r 





The equations for the x and y scale factors S xi and S Y; j_, 
respectively, are derived by assuming that the scale factor 
varies linearly from the bottom to the top for x and from left 
to right for y. For example, we can compute from known data 
a y scale factor on the left, S YL , and a y scale factor on 
the right, S Y p. The x position of these scale factors are known 
data points X* and x R . On the basis of the linearity assumption, 
we can write the general equation: 

S Yi = m Y x i + b Y (17) 

But the two end conditions yield: 


m y - (S YR - S yl )/(x r - x^) 
b Y = S YL “ ®Y *L 


Similarly for S^: 


S Xi ” m x 9 i + b x 


here the two end conditions yield: 


m X ^vt “ I (Yt “ Fi*) 


XT XB 


W 


~ t 


b X “ S XB - m X 


Slope, Slope Length, and Aspect Determination 


(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 


Thus, equations (2) and (3) yield an ,Z£> type 

point for each plate coordinate point. If the output cell size 
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is greater than the input cell of .01" x .01" (208' x 208'), 
then the average of all the Z's in the output cell domain is 
computed and used as the Z value for the cell size. 


For purposes of output and for subsequent entry into the 
data base, a column qnd row number for each point is computed 

from: 

* 

C.N. = (X E ^- Xq) + 1 J 

c e ( (23) 

R.N. - <Y n - Y 0 ) + 1 J 
c e 

where : 

C. = output cell size 

(Xq/Yq) - origin of output area, in meters 


Consider the following figure 8 

1 ° 


J7‘ 


^ h 


X i Y 3 Z 3 

P 6 

x 2 y 3 z 6 

P 9 

X 3 Y 3 Z 9 

P 2 

P 5 

p 8 

W2 

x 2 Y 2 z 5 

x 3 y 2 z 8 

p i 

P 4 

P 7 

x i Y i z i 

X 2 YiZ 4 

w 7 


■IH 


Figure 8.- Neighborhood cell approach. 
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Slope is defined as IaH | / AD, i.e.: 
s 1 ( p 5 P I ) - <z 5 - z x ) / \/x 2 + AY 2 

S 2 ( p 5 p 2 ) = (Z 5 - z 2 ) / (X 2 - Xj) 

S 3 ( p 5 p 3 > » (Z 5 - Z 3 ) / \Ax 2 + aY 2 

S 4 ( p 5 P 4 ) = (Z 5 - Z 4 ) / (Y 2 - Yj) (24) 

s 5 ( p 5 p 6 ) - (Z 5 - z 6 > / (Y 3 - Y 2 ) 

S 6 (P 5 P 7 ) - <Z S - Z 7 ) /v4 X 2 + 4 Y 2 

S 7 ( p 5 p 8 ) » CZ 5 - Zg) / (X 3 - X 2 ) 

Sg( p 5 p g ) - (Z 5 - Z 9 ) /\/aX 2 + AY 2 

The largest (i = 1,2,.., 8 ) is selected as the slope of 
the cell containing P 5 . The Z 5 value and the Z^ corresponding 
to the largest are then compared. Then the aspect is defined 
as the direction from the smallest Z to the largest Z. 


Slope length is given by the following equation: 


S.I.. = v/" 2(C 0 ) 2 + AZ 2 


(25) 


Since the slope and aspect determination are based on a 
neighboring approach, the perimeter cells for the adjacent 
files of a quad map are duplicated for us^ in the algorithm. 
This assures perimeter values for the slope and aspect of a 
1° x 2° area at the tape file interfaces. 
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RESULTS 


S hiprock and Gallup Quad Data 

The digital terrain data for the Shiprock and Gallup Quad, 
located as shown in figure 9, were selected for processing and 
for subsequenr entry in the IDDP. The cross-hatched area in 
this figure represents the area of interest. Input parameters 
such as desired output origin, cell size, quantum levels for 
elevation, and slope were determined. Figure 10 shows a 
typical computer printout for one of the variables (topographic 
elevation) in the data base. The input data for the example 
shown in figure 10 Is given in table I. The input indicates 
241 rows and 354 columns of output, but only 62 rows and 120 
columns are shown for the sake of brevity. Each letter repre- 
sents a cell ICOm x 100m in size. Each cell is identified by 
a row and column number and can be related back to the UTM 
coordinate system by equation (23) . It should be pointed out 
that the (x QJ y 0 ) origin is referenced to the lower left-hand 
cell corner. 

To determine the magnitude of agreement between the actual 
data and data produced from the TOPO program, the elevation data 
for these files were aligned with the corresponding Lands at scene 
data. Specific mountain peaks and/or features for the file 2 
data were then identified from a 1 :62, 500-scale USGS map. A 
northing, easting, and elevation value were read from the map 
for three specific points. The northing and easting of each 


Zone 

CM 



1: File 1 

2 : File 2 

3: File 3 

4 : File 4 


Figure S. Shiprock and Gallup Quad Data, RCC 135 
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TABLE I.- INPUT DATA FOR EXAMPLE SHOWN IN FIGURE 10 


INPUT COLUMNS TO PROCESS - 1 354 

INPUT STARTING ROW ** 1 

CELL SIZE = 100.00 

NUMBER OF ROWS OUTPUT = 241 FROM 1 TO 241 

354 COLUMNS OUTPUT NUMBERED FROM 1 TO 354 


ELEVATION 

LEVELS 

COMPUTED 

IN INCREMENTS OF 

100 . 0 

FROM 6000 

LEVEL 


ST. 

EL. 

END EL 

* 

1 

A 


0.0 


6000 . 0 

2 

B 


6000.0 


6100.0 

3 

C 


6100.0 


6200.0 

4 

D 


6200.0 


6300.0 

5 

E 


6300.0 


6400.0 

6 

F 


6400.0 


6500.0 

7 

G 


6500.0 


6600.0 

8 

H 


6600.0 


6700.0 

9 

I 


6700.0 


6800.0 

10 

J 


6800.0 


6900.0 

11 

K 


6900.0 


7000.0 

12 

L 


7000.0 


7100.0 

13 

M 


7100.0 


7200.0 

14 

N 


7200.0 


7300.0 

15 

0 


7300.0 


7400.0 

16 

P 


7400.0 


7500.0 

17 

Q 


7500.0 


7600.0 

18 

R 


7600.0 


7700.0 

19 

S 


7700.0 


7800.0 

20 

T 


7800.0 


7900.0 

21 

U 


7900.0 


8000.0 

22 

V 


8000.0 


8100.0 

23 

W 


8100.0 


8200.0 

24 

X 


8200.0 


8300.0 

25 

Y 


8300.0 


8400.0 

26 

z v 


8400.0 


99999.0 


point were then used to compute a row and column number and 
the corresponding elevation interval was determined from the 
TOPO output data. This comparison is shown in table II and 
shows close agreement between the actual data and the data 
produced as output from the TOPO programs. 

As a further accuracy evaluation, the four known sheet 
comers for files 3 and 4 of the RCC 135 data and file 3 of 
the RCC 136 data were translated and registered to the UTM 
coordinate system using equations (2) and (3). These UTM 
values were then compared with the known UTM values and the 
differences recorded as shown in table III. 

CONCLUDING REMARKS 

The NASA/ERL has developed a method for processing digital 
topography data that can subsequently be entered in a data base 
to include slope, slope length, elevation, and aspect. It is 
expected that this information, and subsequent second and third 
level interpretive information derived from the original source 
data, can be used by land resource managers. 

The current software programs have been written for an 
input data tape formatted as per the DMA Planar map data file 
format. Another tape format from the DMA is available to users 
whereby the data are referenced to the latitude , longitude 
coordinate system in terms of 3, 1, or .5 arc sec cell sizes.' 
The former data format is used at the ERL since the majority 
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| 

L 

TABLE II.- COMPARISON OF ELEVATION VALUES FROM MAP 
} AND TOPO PROGRAM 


Point Number 

Elevation from 
1:62, 500 map , ft . 

Elevation from 
TOPO Program, ft. 

I 

7000 

7000 - 7100 

2 

8171 

8000 - 8100 

3 

; 

8304 

8100 - 8300 


TABLE III.- UTM DIFFERENCES FOR KNOWN POSITIONS 


























of the Laboratory applications are concerned with the utili- 
zation of map data that are based on the UTM reference grid 
system. 
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APPENDIX A 


NASA/ ERL DIGITAL TERRAIN TOPOGRAPHIC PROGRAMS 

This appendix defines the procedures and technique involved 
in processing digital terrain data tapes produced by the DMATC. 

This processing necessitates the execution of four separate 
programs from which is output four data files (elevation, slope, 
aspect and slope length) which are in a form suitable for entering 
in the Gridded Data Base.' Additional information can be obtained 
in the DMA TOPO Program Documentation manual on file with the 
NASA/ERL. The hardware required is a Varian V70 series mini- 
computer with the following program memory requirements : 

Program 1 - 51100g 

Program 2 - 63600g 

Program 3 - 47122g 

Program 4 - 47373g 

One card reader, a line printer, and two tapes or disk 
files are required. The individual programs are run in the 
sequence as defined in this appendix. 

/ 


AVI 


PROGRAM 1 


TOPREF 


/END JOB . 





/EXEC, T0PREF. 

/ACCESS, CDM0UT, MT , , REFORMATTED OUTPUT TAPE. 



/ COPYF, CDMA IN, ,F1. SKIP FILE(S) TO DATA FILE. 



/ACCESS, CDMA IN, MT, REELN0. N CIC I NPUT TAPE. 



( Front of dtck) 


* Data is contained in the second tape file of each group of 3 files 
for each area. For a four area tape, skip 1, 4, 7, or 10 files to 
be positioned at the desired data file. 


A- 3 








CARD 


NO. 


JOB TOPO REFOREST 
NAME T0PREF 


LEAD CARD SET UP 


PROGRAMMER 


J. Forbes 


PAGE NO.Jl_.OfJ: 


DATE 


5/257W 










PROGRAM 2 


TOPTWO 



(Front ol d»ch) 





CARO 


NO. l 

JOB ..-TP.g.o.i^.SSPW 2 , lead CARD SET UP page N0.i_0FJL 

NAME JPgffjg PROGRAMMER J - Forbes DATE 6/23 / 78 

r r m - '■■■> *» I ■' rnmmm ■■ — i mmmmm » mm*m ■■ i » ■ nn ■ mmm m mm . ^ . ■ ■ ■ - 


(FIELD f CARD 
J 1 D {COLUMNS 

FORMAT 

SYMUGLIC 

name 

IDENTIFICATION 

r 

i i 

\ 1 \ 1-12 

D12.2 

X0 

Easting output origin in meters, minus 

1/2 i 

‘ 

: ! « 

• 


output cell size. 

i 

t 

- ' 

; i j 


1 ; 

1 

i | 

i 2 ! 13-24 

D12.2 

• 

* Y0 

i 

Northing output origin in meters, minus 1/2 

; : ! 


output cell size. * I 

: ; i 

: » : 

i ■ i 

f 1 

• 3 •' 25-36 

D12.2 

XE 

n 

? 

1 Easting output limit, plus 1/2 output 

l» 

cell size, j 

1 “ - * - ■ ■ 1 1 - 1 

i 

i 

i 

i - 

i 

t 

f 

j 

' 4 ‘ 37-48 

D12 .2 

YE 

1 

Northing output limit, plus 1/2 output 

cell size: 

» 1 
1 * 



— 1 


| 

1 

5 • 49-60 

D12.2 

X0RG 


Map sheet corner, easting I 

1 

• ; i i _ J 

6 61-72 

012 . 2 ' 

Y0RG 

t 

J 

Map sheet corner, northinel . 

j 

1 * 






j 






ORIGINAL PAGE IB 
. OF POOR OTTAT.TTV 

-i 

i 

** Field input example ** 


1 



— 

649950. Input as 649950D2 





i 

1 

3985597.81 Input as 398559781D0 



. 


■ 

1 


> 




! 

1 

1 





1 


| 

COMMf- ht^ ^ Map sheet corner 

= 

SW corner for east file and SE corner for west 


fi le of qua d map. 











CARD 

NO. 2 

JOB TOPO - PROGRAM 2 

NAMF T0PIW0 


LEAD CARD SET UP PAGE NO._2_CF_6 

PROGRAMMER J - Forbes DATF 5/25/73 


FIELD CAHD SYMBOLIC 

I. D. COLUMNS ' OR * AT NAME 


IDENTIFICATION 


15 ISTA.RT Record number on input file to start processing 


Deafult 


_ 2 

6-15 

F10.0 

CELL 

Output ^11 size in meters. 


3 

16-25 

F10.0 

CB 

Column bias in meters. 


26-35 | 

J 

E10.0 

1 

RB 

Row bias in meters. 


COMMENTS 


























CARD 

NO. 

j 0 B TOPO-PROGRAM 

NAME - TOPTHO 


CARD' 

COLUMNS 

FORMAT 

SYMBOLIC 

NAME 

• 1-20 

D20.10 

DELTA 

21-40 

D20.10 

THETA 

41 

11 

LQ 






LEAD CARD SET UP ' PAGE NO._3_OF.6_ 
PROGRAMMER Jt Forbes ... DATE._jZ 25/78 

IDENTIFICATION 


=1 for processing left quadrant, otherwise blank. 


■^The DELTA angle is determined from: 
tan A = J E 2 - E-^ | /(N2 - %) where the 


subscript 2 refers to the top sheet coordinates 


and the subscript 1 refers to the bottom sheet 


coordinates of the Quad map centerline. For an 


input data file relative to the Zone CM, the 


"sign" of the DELTA angle is- as follows 


—Negative: East of Zone CM, left half of Quad 


—Positive: East of Zone CM, right half of Quad 


—Positive: West of Zone CM, left half of Quad 


—Negative: West of Zone CM,- right half of Quad 



“THETA is normally zero, but should be 


verified by checking the X sheet corners 


of the central meridian of .the Quad map. 










CARD 

NO. k 

JOB TOPO - PROGRAM <2 

NAME 


TOPTWO 


LEAD CARD SET 

PROGRAMMER . 


UP PAGE NO.. 

J. Forbes DATE 5/28/78 










CARD 
NO. _ 

joa - 

NAME 


TQPO - PROGRAM 2 
TOPTW0 


LEAD CARD SET UP 

PROGRAMMER J - I°-£. b ^ s .. 


PAGE NO.JL.OF_L 
DATE 5/28/78 






















CARD 

NO. 6 

JOB -TOPO gROgSAMj LEAD CARD SET UP PAGE NO._6_OF_jL. 

NAME .12 PROGRAMMFR J - Forbes nATF 5/28/78 



I 


A- 12 . 






















PROGRAM 3 


/ 


- TOPSRT 


( 


/ENDJOB 




DATA CARD 


T A 


/EXEC, T0PSRT. 




CARL) 

NO. i 

;(■!! Jl li(X L ! <AlL? LEAD CARD SET UP PAGL n o_Lof.X 

name PROGRAMMER DA TE M^/l IS— 


ilRlof CARO ....... 

i o 1COIUUNU * 0MiA ’ 

liYMUOUC 

wam£ 

IDENTIFICATION ! 

1 ' 1-5 15 

NREC 

Number of records on input file 1 CDMiflOT ' . 

i 

* i 


! 

2 i 6-10 ! 15 

* 1 

— — - — n 

IR0W 

_ 

, 

Starting row number. 

- - ir - ■ , ■ .. ■ — ^ 

; I 

i ! ; 

; ; < 

i 

1 

, 3 11-15 15 • 

LROW 

a 

Last row number. 



, , , i 

4 I 16-20 15 ! NT ROWS ! Total rows (calculated in program 2). 



5 21-25 15 i NCOL ; Number of columns. 



This input is determined from printed output f rom previo us F!i(.x;_R/d;I-2 


A- 15 






PROGRAM 4 


TOPODB 


DECK SETUP FOR PROGRAM TQPO ...» PROGRAM 4 PAGE NO. 

(Bock of dock) 


L 


f 


/ENDJOB 


( 


(. 


( 


f. 




DATA. CARDS 


/EXEC, . T0P0DB . 


( 


(. 


! /ACCESS, LNDATA, DP. 2 


/ACCESS, ASDATA, DD. 2 




/ACCESS, SLDATA , DP.' 


/ACCESS, ELDATA, DD. ' 


• / 


. ( 

( 

| /CREATE, SLDATA, 1080, 


/CREA TE, LNDATA, 1080, NRECS/S. SLO PE LENGTH DATA. 2 
/CREATE, ASDATA, 1080, NRECS/S. ASPECT DATA. 2 

NRECS/S. SLOPE DATA. 2 


/CREATE, ELDATA, 1080, NRECS/S. ELEVATION DATA. 


(l 


a 


/ACCESS, SRTTAP, MT, L9XXXX, OUTPUT FROM PGM-3. 


/ JOB, STANDARD JOB CARD. 


(Front of d»ck) 

* Number of records output is determined by number of rows processed 
(see Lead Cards). 

1 Omit this file if I0PT = 2 (lead card #1) ' 

2 Omit these files if I0PT = 1 

A-17 


«■ »» «>>. mt . - n — m . n l i 


CARO 

NO. 1 

J0B TOPO - PROGRAM 4 

NAME TOPO DB __ - 


LEAD CARD SET UP PAGE NO.-l_OF.-J 

PROGRAMMER J . Forbes DATE 5/28/78 

























OPO - PROGRAM 4 
TOPODB 


LEAD CARD SET UP PAGE NO.JL.OE_Z 

PROGRAMMER J., Forbes DATE 5/28/78 


field caro symbolic 

t. o. columns ruh NAME 


identification 


5 15 IROW First row number for output file (input to 


database). 


2 

6-10 

1 

15 , 

i 

| LRCW 

Ending row number fbr output. 


11-15 j -5 | ICBIA.S | Column bias for shifting output. 





















CARD 
NO. _ 


TOPO - PROGRAM 4 


NAME - T 0PQD 1. . 

Field j " card 
I. D. COLUMNS 


! FORMAT 


SYMBOLIC 

NAME 


LEAD CARD SET UP PAGE no.JL.of 7 

PROGRAMMER — J ' Forbes DATE _ 5 / 28 J Z*t 

IDENTIFICATION 


i 1 1-8 F8.0 STELV 


Starting elevation level. 








JOB gQPQ PROGRAM. 4 LEAD CARD SET UP PAGE NO..*_-Or_JL 

MAMF T0PQDB ; PROGRAMMER DATE 5/2S/ ?L 


( FIELD 
i I D. 


CARD 

COLUMNS 


format 


SYMBOLIC 

name 


IDENTIFICATION 


1- 5 


15 


NSL 


Number of slope levels Co read (following), 


MAX - 20. 



i 



Omit this card If I0PT ■ 1 


A- 2 1 


COMMENTS 


CARO 

NO. - 5 

JOB T0PQ - PROGRAM 4 

NAME TQPQDB - 


Field 
». o. 


LEAD CARD SET UP page no. JL_of_Z 

PROGRAMMER DATE 



ff'ffifflfflimIBPBi 

IDENTIFICATION 1 

1- 8 

F8.0 

XSSL (I) 

Starting slope level (1). ! 



F8.0 XSSL(3) 


25-32 F8.0 XSSL(4) 


33-40 ! F8.0 XSSL (5) 


6 41-48 F8.0 XSSL (6) 


S 7 149-56 ! F8„0 ! XSSL (7) 

> I ^~1 | 

f 8 57-64 F8.0 XSSL (8) 

? 1 


! 9 '65-72 I F8 t 0 XSSL (9) 


Starting slope level (2). 


Starting slope level (3). 


Starting slope level (4). 


Starting slope level (3)'. 


Starting slope level (6). 


Starting slope level (7)« 
Starting slope level (8). 
Starting slope level (9). 




















4 LEAD CARD SET UP PAGE NO. 

PROGRAMMER For be s DATE 


IDENTIFICATION 


IASP(l) Aspect number code for undetermined. 


IASP(2) Aspect number code for NORTH. 


12 IASP(3) Aspect number code for NE. 


12 IASP(4) j Aspect number code for E 


12 IASP(5) jAspect number code for 'SE 


IASP(6) jAspect number code for S 


12 IASP(7) {Aspect number code for SW 


12 IASP(8) jAspect number code for W 


12 IASP(9) jAspect number code for NW 


COMMENTS Omit th is card if I0PT =_1 



































CARO 

NO. Z 

JOB TOPO - PROGRAM 4 

NAME 


TOPODB 


LEAD CARD SET UP 

PROGRAMMER J - 


PAGE NO.. 
DATE 


5/28/78 

















g»«Jg 


Mr* 


APPENDIX B 


jj . 

» * 


PROGRAM FLOW CHARTS 



* 



B-l 



4 




Program 1 - TOPREF 

" ■ ■ ■ f»— ' 1 ' - ' 













Program 3 (continued) 





PROGRAM 4 - TOPODB 



(continued) 

B-6 




Program 4 (continued) 




Read disc files 
print out data. 




4 


Print | 
elevation I 
frequency j 
distribution! 




PRINTER 

MkPS 




AT 

S 

' APPENDIX C 

;; | 

PROGRAM LISTINGS 



PAGE T 01/04/78 47103000 VORTXII rtf 'IV(G) A148AJ 0921 HOURS 


I 


NAMF TOPREF 

3 

1 

2 


TITLE M48 A l 

3 

2 

3 


DIMENSION INC 1 SR40) ,TA (2005) ,IFN (3) 

3 

3 

4 

C 

„ • . ... 

3 

4 

5 


INTEGER OFNC37 

3 

5 

6 

c 


3 

6 

7 


DATA IFN ,OFV/ • CD • , ' MA * TN ' , f CD 1 , ' MO* , f UT 1 / * 

3 

7 

8 


DATA IBLK/» '/ 

3 

0 

9 

1000 

FORMATC2I6, 315,9X316, ' .,f',l5,» )..',3X6) 

3 

9 

10 

1001 

PORMATC ' FLAG REC X’,4X*T1 YN'fllX'ZJ Z2 Z3», 

3 

10 

11 


X 7X*N%7X’ZN-2 ZN-I ZN'7 

3 

11 

12 

1005 

FORMAT(lHl) 

3 

12 

13 

1006 

FORMAT (// , ' MIMfMAX ELEVATIONS, *,2110) 

3 

1 3 

14 


CALL CSINTT 

3 

14 

15 


CALL PAGES 

3 

15 

16 


WRITEC5, 1Q0I) 

3 

16 

17 


MIN = 32000 

3 

17 

18 


MAX * -9999 

3 

18 

19 


KR - 0 — ’ 

3 

19 

20 


IA (2 ) * 0 

3 

20 

21 


I Flu * 0 

3 

21 

22 


TA ( 1 ) * irL 

3 

22 

23 

90 

CONTINUE 

3 

23 

24 


DO 95 I • 1,15840 

3 

24 

25 

95 

IN ( I ) * 0 ” 

3 

25 

26 


CALL VSRrdFN, 15840, IN, L) 

3 

26 

27 


IF fL ,GT, 0) GO TO 98 

3 

27 

28 


IFCL , EO. -3) GO TO 900 

3 

28 

29 

98 

CONTINUE 

3 

29 

30 


LC * IN(9) 

3 

30 

31 


NWDS - L • 

3 

31 

32 


KR s KR + 1 

3 

32 

33 


IF ( KR .EO. 1) GO Tft 115 

3 

33 

34 


K * 5 

3 

34 

35 


GO TO 120 

3 

35 

36 

115 

K * 46 

3 

36 

37 

120 

CONTINUF - - ~ 

3 

37 

38 

C 


3 

3R 

39 

130 

CONTINUE 

3 

39 

40 


IF ( IN (K ) .EO. -17 GO TO 90" * ‘ 

3 

40 

41 


DO 155 *L * 6,2005 

3 

41 

42 

155 

IA ( L ) = 0 

3 

42 

43 


I = INCK + 4) " " * " ” ' * 

3 

43 

44 


T A C 2 3 = I A ( 2 ) + 1 

3 

44 

45 


TA ( 33 = INCK) 

3 

45 

46 


1 :3 C 4 ) « TNCK+1) 

3 

46 

47 


X A C 5 ) = I 

3 

47 

48 


IVN * I A ( 4 ) ♦ I - l 

3 

48 

49 


NL * K*5 - - 

3 

49 

50 


NU * K+I+4 

3 

50 

51 


N - 5 

3 

51 

52 


DO 160 LsNL.NU - - 

3 

52 

53 


N = N * .1 

3 

53 

54 


IACN) * TN(L) 

3 

54 

55 


IF fTA c N ) .GT . WATT~»ntT*I ATrT ) ~ “ " 

3 

55 

56 


IF ilA(N).LT.MIN) MINsIACN) 

3 

56 

57 

160 

CONTINUF 

3 

57 

58 


NRITEC5 , 1 0007 tTAtb) ,L»1 /tTrlTW. CIA CL7 ,L*6 , 8 ) , 1 2IA tN-27 , IACR-1 3 

3 

58 

59 


X , I A (N) 

3 

59 


C-2 ORIGINAL PAGE IS 

OE POOR QUALITY 


p agf 


2 0t/04/7H 47193000 VORTXII 


r TM I V C G > A14IA1 


0923 HOURS 


60 

165 

CAM- V$ WR (OPN , 2005, IA »L) 



3 

oO 

61 


IF Cl, . F.O. -2) GO TO 165 



3 

61 

62 


K « »U(IN{K-2)/2! 

"* * ' ’ **‘ i ’ 


3 

62 

63 


IF(K , GT. NWOS) GO TO 90 



3 

63 

64 


GO TO 130 


. 

3 

64 

65 

C 

' 

m, *■ .. ~ 


3 

65 

66 

500 

CONTINUE 



3 

66 

67 


WRITE ( 5,1 006 ) HIM, MAX 



3 

67 

66 


WRTTEC5, 1005) 



3 

68 

69 


CALL VAEF(OFN) 



3 

69 

70 


CALL VIEF(OFN) 



3 

70 

71 


CALL VIRF-vlFM) 



3 

71 

72 


CALL V«RE(OPM) 



3 

72 

73 


STOP 



3 

71 

74 


END 

* — - ' t* <•—* 


3 

74 






C-3 


PAGE 

1 

2 

3 

4 

5 

6 

7 

8 
4 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
2 3 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 
34 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 
-55“ 

56 

57 

55 
. 59 

i 

i 


3 01/04/78 47183000 VORTXTI FTN IVCG) B148A1 0924 HOURS 

NAME TOPTWO 
TITLF ft 1 48 A 1 
INTEGER f)FN 

DOUBLE PRECISION XO, YO,XE,YE,XORG, YORG, XI ,Y 
DOUBLE PRECISION YM , B Y , XM , BX , CON , DBL l , DBL2 
DOUBLE PRECISION XC, YC,XP,YP,SX,SY,XPP,YPP 
DOUBLE PRECISION SYL,SYR,SXB,SXT,XOR,YOR,XCELL ' " 

DOUBLF PRECISION DELTA , THETA , YC1 

DIMENSION TFN!3),0FN! 3),Y!2000),YA!2000),IA(2005),IB!7200) 
EQUIVALENCE !XFLAG,TA(1)), (IRr IAC237 VtLXf IA(3) ) 

X ,(LY,IA!4)),!N,IA(5)) 

DATA CON/O. 017453292D0/ 

DATA IFN , OFN/ ' CD ' , »MA • , ’ IN' , rCD HO* , »UT'/ 

1000 FORMAT! • OUTPUT UTM ORIGIN ',2015.9/ 

X » UPPER RIGHT LIMITS' , 2015.9,/ 

X • DMAATC UTM ORIGIN »,2D15.9/) 

1001 FORMATI6D12.2) 

1002 FORMAT! 15 , 3F1 0 , 0 ) 

1003 FORMAT!/' RUN COMPLETE - WROTE' ,15,^ RECORDS*/, * 

X' FIRST ROW = ' ,15, • 

X* FIRST COLUMN «• ,15, ' 

1005 FORMAT!/' FIRST X ! ,I6,' 

1006 FORMAT! 415, Srt5. 2* 

COL 


LAST ROW * 

LAST COLUMN 
FTRST T'.Tin 


' ,15/, 
' ,15/) 


OUT 


J' , 12X'XV, l IX' Y(l) ' ,tlX't(l)',llX'Y(J)', 


1007 FORMAT ! 3 X ' X M 
X UX'Z!J)»/7 

1010 FORMAT!2D20.10,I1) 

1011 FURMAT!2F10.0) 

1012 FORMAT C4D20 .10} • 

1013 FORNAT(4F10.0) 

1100 F0RMATC1X, 'CALCULATED TOTAL NUMBER OF ROWS OUTPUT s;,15,/, 

♦ IX, 'CALCTEATED iWMtElTOF-COLtOrNS-OOTPtrr »^TS ,T, 

* IX, 'ROW TO START PROCESSING ',15,/, 

♦ IX, 'OUTPUT CELL SIZE ',F6.2,' METERS',/, 

* ix, 'column bias* ;tio.27/7 

4 IX, 'ROW BIAS ' ,F10,2, ' METERS') 


DELTA * ',020.7,/, 

THETA* ♦ ,D20;7-,/, 

QUAD FLAG * ' ,Il,/, 


1102 FORMAT!' 

* 

* 

* 

4 

1103 FORMAT! 

* 

9 
•* 

• 

OLL -CBTHTT 

CALL PAGES 

READ! 4 ,1001)XO,YQ,XE,YE, XORG , YORG 

RE A P! 4, 10 02) I 3 T A RT , CELL, C », R » 

READ! 4, 1010) DELTA, THETA, be 
REAp(4,!0il) XOAT,YORI 
RE A D ! 4, 101 2) S T L, 


XORZ * 
TORT *r 

SYL * 
SYR * 
SXH 3* 
SXT = 
XL s 
XR * 
YB « 
YT ■ 


' , E 1 2 . 2 , / , 

nrmi) — 

',020.10,/, 

' , D20 . 10 , / , 

*, 020710“,/; 

' , 1 ) 20 . 10 ,/, 

» , F12 . 2 , / , 
FTT,2,rr 

•,m.2,/, 

»,ria,2) 




S TR ,S X B,S XT 


READ !4,1013) XL, XR,YB,YT 
If (I ST ART.EQ.O) ISTART«1 

WCOL ■ (XE*XO)/CELL 


9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

• 

5 

9 

9 


1 
2 

3 

4 

5 

6 
7 

e 
s 

10 
n 
12 

13 

14 

15 

16 
17 

ie 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 
3i 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
41 
19 

50 

51 
9 52 

• S3 

• 54 
-9~ 55 

9 56 

9 57 

“5 St 

• 59 


C-4 


f 


01/04/78 47133000 VORTXII FTN IVC<3) B148A1 


0924 HOURS 


60 


WRITE! 5, 1000 ) XO,YO,X£,YE,XORG,Yr»RG 

9 

60 

61 


WRITE! 5, 1100) NRQW ,NCOL » ISTART , CELL ,CB , RB 

9 

61 

62 


WRITE! 5, 1102) DELTA,THrrir7IiQr,inirr;TORI “ * 

9 

62 

63 


WRITE! 5,1103) SYL.SYR ,8XB;8XT,XL,XR,YB,YT 

9 

63 

64 

C 

* 

9 

64 

65 


LROW = NPOW 

9 

65 

66 


KROW = 0 

9 

66 

67 


XCELL 9 CELL 

9 

67 

66 


THETA ■ V't&?4*C3N 

9 

68 

69 


DELTA • LELTA«CON 

9 

69 

70 


OBL1 • CB 

9 

70 

71 


XO * XO+DBLt * 

9 

71 

72 


XE * XE+DBLt 

9 

72 

73 


DBL1 * RB 

9 

73 

74 


YO * YO+DBLl 

9 

74 

75 


YE « YEtDBLl 

9 

75 

76 


DBL1 9 XR 

9 

76 

77 


DBL2 * XL 

9 

77 

78 


YM = (SYP-SYD/CDBL1-DBL2) 

9 

78 

79 


BY = SYL-YM4DBL2 

9 

7 V 

80 


DBL1 « YT ' 

9 

80 

31 


BBL2 « YB 

9 

81 

82 


XH « (SXT-SXB)/(DBL1-0BL2) 

9 

82 

83 


BX * SXB-XH4DBL2 " ^ * 

9 

83 

84 


ICN 9 0 

9 

84 

85 


IVL = 0 

9 

85 

36 


NRO * 0 ’ "• 

9 

86 

87 


IXC « 1 

9 

87 

81 


A ■ 0.0 

9 

88 

89 


Ycr « o . our - - • 

9 

69 

90 


I YP = 0 

9 

90 

91 


IYL = 1 

9 

91 

92 


NLI* "» 65 

9 

92 

93 

• 

DBL1 ■ D5INCDILTA) 

9 

91 

94 


DBL2 > DCOSC DELTA) 

9 

94 

95 

200 

CONTINUE 

9 

95 

96 


CALL VSRDUFM, 2005, IA, 1ST) 

9 

96 

97 


IF! 1ST .GT. 0) GO TO 205 

9 

97 

91 


IF CIST , CO. -2) GtrTD 2*0 T ‘ ~ 

9 

98 

99 

205 

CONTINUE 

9 

99 

100 


IFCLX .LT. 1) GO TO 200 

9 

too 

101 


ICN -B TCN ♦ 1~ * ' ■ “ ■ ~ 

9 

101 

102 


IFCICN .LT. ISTART) GO TO 200 

9 

102 

103 


IF t ICN .GT. ISTART) GO TO 220 

9 

103 

104 


XC « LX " “• ■*" ■ ■ 

9 

104 

109 


YC b LY 

9 

105 

10* 


WRITE!*, 1005) LX,LY 

t 

106 

TOT 

- * “CREtr pxstr* . - ----- - - 

■ ‘ - 9 

107 

108 

220 

CONTINUE 

9 

108 

109 


XP * LX 

9 

109 


-TTTTC ' 

XOR ■ XORZ 
TOR " TORI 

Xp - B - X paq t O B — 

YP * YP-YOR 

IF C LO ) 230,225,2 30 

tlf |FT r (XP*nin *~CYP*C 0*1 
ipp • (xp«ma>”(YP*DBU> 


9 110 
9 ill 
9 112 
~9 111 
9 114 
9 US 

-f 11 * 

9 117 
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118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 


230 

235 


240 


252 

255 


GO TO 235 

YPP « C-XP»DBLn - TTFnPftn « P »bn/DtL 3r ♦CTP/PILl) 
XPP > (XP*0BL2)4(YP*DIL1) 

CONTINUE 

SX s XN«YPP+BX * • • 

SY * YM*XPP*BY 
XI * XPP4SX+X0RG 

IF CXl.LT.XO) GO TO 200 ■ 

IF CX1.GT.X8) GO TO 2«0 
Y f{ 1 > « YPP*SY4Y0RG 

YACt) a XA (S ) • " - - " 

N1 * 7 

DO 240 X a 2,N 

YCI) ■ Ytr-UfSY 

YACI) a XACIfl) 

N1 ■ N1 ♦ 1 

continue - * — - 

IVL * IVL+l 

J a 0 

DOTS 2 I * T 7 » * 

IFCYCI) .LT. YO) GO TO 212 
IF(YCI) .GT. YE) GO TO 2S9 

J x j * 1 “ — 7~ 

YCJ) a Y C X 3 
YA(J) a YACn 

CONTINUE • 

CONTINUE 

IX a (Xl-XOJ/CELL ♦ 1. 


if tix. LETtrn GrrTir-m 

IF (IX.GT.NCOL) GO TO 200 
IFtNLIN .LT. 55) GO TO 261 

CALLTAGEr * ; : 

NRITEf 5, 1007) 

NLIN a 0 


~26T CONY TRUE - 


riix-nc 
YJJ • t(J 



N ti TN - a -Hti r *~t 

WRITEC5/ 1006) IR , I VL , IX , J ,XX , Y1 1 , Y12 , YJl , YJ2 
IFCIVL .EQ. 1) ICOL a IX 

DO 270 X'ni 17U : 

IV • CYCI) - TO)/XCC|ib l, 

SfCTT *LE. 0) GO TO 270 

TFTTT ~ . tT~NRO vn GO TO 200 

IFCIYL .EQ. 1) IRON a IY 
IFdYL.EQ. 1 .OR. XU.EQ.2) GO TO 262 

i r m . EE . I TP) GO yo 2 45 • 

CONTINUE 

A ♦ AB8(YA(f )) ♦ .1 


TT L ■ 3 
IXP * IX 
XYP « IY 
C l ■ TC I 
TO 270 


♦T. PO 
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265 YCS = YC1 

TF (IYP.GT.KPOH) KROW«IYP 

IB(TXC+2) = 4 /YCS * 

IB(IXC.tl) s IYP 
IB(IXC) = TXP 

4 * 0.0 - ' 

yci = o.ono 

TXC = 1XC + 3 

IF(IXC .LE. 7200) GO TO 262 
CALL V8WR(OFN,7200,IB,tST> 

IXC = 1 

IF (1ST ,EQ. -4) GO TO 2*0 ' 

NRO = NRO * 1 
GO TO 262 

270 CONTINUE — — 

IYL * 2 

IF (IYP.GE.LROW) G(T TO 205 

NXTP * IYP 

DO 275 K = NXTP,LROW 

IBCIXC) s TXP “ ~ • 

IB(IXC+1 ) s K 
IB( IXC+2 ) * 0 

IXC s IXC+3 - — • 

IF (IXC. LE. 7200) GO TO 275 
CALL V«WR(OFN,7200,IB,I«T) 

IXC * 1 > : * - 

IF (IST.E0.-4) GO TO 210 
NRO = NRO+1 

275 CONTINMr " ; - - 

GO TO 200 
280 CONTINUE 

IF (TX.GT.7rC0t,r IX * -TT«T T — - 

LCOL * IX 

IF (IXC .EQ. 1) GO TO 290 

00 285 I * -TXC, 7200 ” — 

285 IB(I) = 0 

CALL V8WR(ijFN,7200,IB # IST) 

*"■ W’*‘inro 7 i 1 — ■ — — - — 

290 CONTINUE* 

CALL VIEFCOFN) 

ITRTTFrtSTTtnm nwi;tr ow . FRO M , I CO L , LCOL — 

STOP 

END 
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NAME TOPSRT 
TTTfiE Cl 48 A I 
INTEGER OD AT , OFN , OST 

DIMENSION IOATC72001,0DATn3X>0TJimm;onrm 
DIMENSION AR(108C),IFAR(J) 

DATA IFN , OFN/ 'CD ! »HQ* , »UT' , 'C0» , * NS' , »RT'/ 

DATA IFAR/*5R' , 'TT 1 , ' AP r /“ ~ ' ' ~ 

1000 FORMATC5I5) 

1001 FORMAT ( / ' NO. RCDS. ON O/P DISC FILE' f l5/J 

1002 FORMAT(20X , 16) 

1003 FORMAT!/* NO. RCDS. ON O/P TAPE FILE', IS/) 

1004 FORMAT! 'i EXIT - SEARCHING FOR ROW VALUE »,|4,' ON INPUT RECORD 

* «,I5,/) 

NUMBER OF RECORDS OUTPUT, IN ERROR. NO TAPE OUT.') 
INPUT PARAMFTERS -' ,/, 

NRECS r ,X5 , If ‘ 

FIRST ROW ',15,/, 

LAST ROW ' ,X5,/, 

OF TOTSX ROWS" *71577, 

NCOLS ',15,/J 

1007 FORMAT! 1 OX, ' MOD 10 ROWS WRITTEN TO SRTTAP -•) 

1008 FORMAT C ' 0 NORMAL EXIT, PART If «C*RCHIWG~TOR~!a¥ VALUE t',I4,') 

* ON INPUT RECORD ',X4,/> 

1009 FORMAT!' WRITE CDMSRT REC ',I4,« LAST IN P UT R OW » ',!•) 

CALL CSINIT 

READ (4,1000) NREC, IROW ,LROW,NTROWS,NCO L 

ISKP^* ! NTRQWS-LR0W^4^*^^^*^^^^*^^^ 

NRO » 0 

TOUTT *• 0 * * * 

NROW s LROW-IROW+1 
NSK = NROW * 3 

KK « LROW ♦ 1 5 — ' — 

H * 1 


♦',14,' N 
1100 FORMAT ( ' 0 
1006 FORMAT C ' 1 

4 ' 

« t 

4 * 

4 i 

* I 


DO 300 K * IROW,LROW 


K = K 

|jk ■ - j ■ — 

N « (NROt4lROW-K)43 - 1 
IF ( N ,GT. 7200) GO TQ. 95 

WRRD = 0 

CALL VSFCCIFN.0,1) 

CALL V$RDCIFN, 7200, IDAT, 1ST) 

VRRO * NRRO + t 

IF(IDAT!N) . NE c KK) GO TO •) 
70 CONTINUE 

■* ourTtwr “S’ TOFT CN-TI ' 

ODAT(M+i ) = IDAT (N) 

ODAT ( M+2 ) s IDATIN+1) 

-tr*-**!-- - 

IF(N ,LT. 7200) GO TX) 90 
10 CALL V9WR(0PN,7200#ODAT,0fT) 

IF COST - ; g q. -4) G O T O 3 50 

IF (OST .EQ. -2) GO TO 80 
NRO « NRO ♦ 1 

TRITE (5, T0(T9TWRG,R 

IF CIOUIT.EQ.I) GO TO UO 
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16 
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21 
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28 
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29 
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31 
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34 
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36 
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37 
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38 
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40 
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41 
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42 

j 

43 

2 

44 

2 

45 

2 

44 

2 

47 

2 

48 

1 

49 

2 

50 

2 

51 

TT 

52 

2 

53 

2 

54 

- r 2T 

1! 

2 

56 

2 

57 


rn 
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60 

M = 1 

2 

60 

61 

DO 85 I = 1,7200 

2 

61 

62 

85 ODATCI) * 0 * 

2 

62 

63 

90 CONTINUE 

2 

63 

64 

NsNfNSKTISKP 

2 

64 

65 

IF CN. IT. 7200) GO TO 91 

2 

65 

66 

IF CN.GT.7202) GO TO 95 

2 

66 

67 

N s N-3 

2 

67 

66 

91 CONTINUE 

2 

68 

69 

IFCIDATCN) ,E0. KK) GO TO 70 

2 

69 

70 

93 CONTINUE 

2 

70 

71 

LL « N - 30 

2 

71 

72 

IFCLt. .LT„ 2) LL * 2 

2 

72 

73 

LU * N ♦ 30 

2 

73 

74 

IFCLU .GT. 7199) LU « 7199 * " ’ 

2 

74 

75 

DO 92 I ■ LL ,LU, 3 

2 

75 

76 

IFCIDATCI) .NE. KK) GO TO 92 

2 

76 

77 

IDIF = I-N - 

2 

77 

78 

NSK s MSKFIDTF 

2 

70 

79 

N = I 

2 

79 

80 

GO TO 70 .. . . 

2 

80 

81 

92 CONTINUE 

2 

81 

82 

IFCLU .EQ. 7199) GO TO 94 

2 

82 

83 

IF CKK.NE.IROW) GOTO 9fr V' 

2 

83 

84 

WRITEC5, 1008 ) KK , NRRO 

2 

84 

85 

97 IOUIT = 1 

2 

85 

86 

IF CM-1) 350,350,80 “**■ 

2 

86 

87 

96 CONTINUE 

2 

87 

88 

If CNRRD.EO.NREC) GO TO 295 

2 

88 

89 

WRTTET5 * 1 0047 “Kit, NRRO PT " ” - 

2 

89 

90 

GO TO 97 

2 

90 

91 

C 

2 

91 

92 

94 CONTINUE ' ‘ : ~ “ ~ - - 

2 

92 

93 

N * 7702 

2 

93 

94 

95 CONTINUE 

2 

94 

95 

NRRO - NRRD F T ” ' 

2 

95 

96 

IFCNRRD .GT. NREC) GO TO 295 

2 

96 

97 

N * N-7200 

2 

97 

98 

CALL V*RD(IFN7T2ff07TO*T,TlT7 

2 

91 

99 

GO TO 9i 

2 

99 

100 

295 CONTINUE 

2 

100 

101 

300 CONTINUE - 

2 

101 

102 

350 CONTINUE 

2 

102 

103 

CALL VIEFCOFN) 

2 

103 

104 

WRITEC5", rOOT) HRO ~ " ~ 

2 

104 

105 

NRR5 • NREC-5 

2 

105 

106 

IE CNR0.GE.NRM5) GO TO 37B 

2 

100 

107 

rWWTTE CT",TTO07 ' 

2 

107 

106 

GO TO 800 

2 

108 

109 

C 

2 

109 

110 

3T5 WMTTET5 , 10077 “ 

- , 

110 

111 

CALL VIFCCQFN ,0,1) - 

2 

111 

112 

NRT * 0 

2 

112 

113 

IC ** 0 

2 

113 

114 

400 CONTINUE 

2 

114 

115 

CALL VSRDCOFN,7200,OOAT,OST) 

2 

115 

116 

JPCOSI.CO,*) .Oft. GO TO 400 

“2 

116 

117 

00 500 I • 1,7200,3 

2 

117 
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118 


IC = TC 4 1 

119 


IFfIC .GT. 1) GO TO 410 

120 

403 

CONTINUE 

121 


IY * ODAT(Itl) 

122 


TF( I Y .EQ. 0) GO TO 500 

123 

405 

CONTINUE 

124 


IX * QDATCT ) 

125 


SUM = 0, 

126 


XN * 0. 

127 

410 

CONTINUE 

128 


IF (ODAT ( 1+1 ) .EQ. 01 GO TO 500 

129 


TF( I Y . NE. QDATITfl 1 ) GO TO 450 

130 


IF( IX .NE. OOATCD) GO TO 420 

131 


XN * XN+1.0 

132 


SUM « SUM + ODATCI+2) 

133 


GO TO 500 

134 

420 

CONTINUE 

135 


APCIX+1) = SUM/XV 

136 


GO TO 405 

137 

450 

CONTINUE 

138 


AR(IX4I1 « SUM/XN 

139 


ARID * IY 

140 


MO * MODCIY ,10) 

3 41 


IF tM0.FQ.0l WRITE! 5 , t 0027 IT 

142 


CALL V$WR(TPAR, 2160, AP. 1ST) 

143 


NRT * NRT + l 

144 


DO 480 J * l,NCOL 

145 


ARCJ) * 0.0 

146 

480 

CONTINUE 

147 


GO TO 403 " 

148 

500 

CONTINUE 

149 


TC = 1 

150 


GO TO 400 • “ 

151 

o 

o 

* 

CONTINUE 

152 


ARCIX+1) * SUM/XN 

153 


ARCH = IT 

154 


WRITE (5,1 002 } IY 

155 


CALL VSWRdFAR, 2160, AR, 1ST) 

154 


NRT * NRT ♦ t ‘ 

157 


WRITE*(5, 1003) NRT 

151 


CALL ViErUPAR) 

159 


CALL “VSEFITFXRT 

160 

800 

CONTINUE 

161 


STOP 

162 


END 
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NAME TOPDDB 
TITLF 014841 
JNTEGFB EFN ,SFN , 4FN 

DIMENSION ARC 1080) ,IXELV( 1080), IXSLPf 1080), IXASP( 1080) 

DIMENSION ELEC(1000,3) ,AP(0),Dirr(8),DIST(8),ICH(2O),IASP(9) 
DIMENSION PARM (2 ) , ELE ( 2) , SLQC2) , ASP( 2) , IFNC-3) , EFN ( 3 ) , SFn( 3) , AFN( 3 ) 
DIMENSION X SSL (20) , XFSM20) , XSCL(254) , XCEL(2S4) 

DIMFNSION TXLFN(1080) ,LEN(3),SLEM(2) 

DIMENSION FLF(254) 

DATA ELF/25440.0/ 

DATA LFN/'LM', 'DA' , 'TA'/ 

DATA SLEN/'S.LE* , 'NOTH' / 

DAT* FLE*SLO,ASP/'FLFV* f * . * , 'SLOP 
DATA IFN ,EFN , SFN , AFN/ 1 SR 1 , 'TT' 

X , 'AS' , 'DA' , »TA'/ 

FORMATUH ,2A4,' DATA COLUMNS' ,TS, * 


E'CASPE*. 'CT'/ 
AF»,'EL','DA',*TA','SL» 


'DA* , 'TA' 


TO 


,15) 


10) •) 


FORMATC (DIVIDED BY 
FORMAT(1XI4,1X120A1) 

FriRMATnXI4,lX120in 
FORMA T( 15, IX, 3014) 

FORMAT(3I5,F10.2,4X,T 1 ) 

F0RHATOI5) 

FORMATdOFS.O) 

FORMATUOI2) 

FORMAT ( ' 1 INPUT COLUMNS TO PROCESS ** » ,217 , /, 

' INPUT STARTING ROW * ',15,/, 

’ CELL 5IZE * ' ,F8 ,2 , / , 

' PROCESSING OPTTOr « ^711,77, 

• NUMBER OF ROWS OUTPUT ■ ',15,' FROM 
15,' COLUMNS OUTPUT NUMBERED FROM ',15,' 

1199 FORMATC OELEVATTON LEVELS COMPUTED IT! INCREMENTS 

* ' FROM • , F8 . 1 , ' TO ' ,F8 , 1 , / , 

* ' LEVEL' ,16X, 'ST. ELV. ' , 15X , 'END ELV. ' , /) 

1200 FORMATC NUMBER OF QUANTUM LEVELS TOR SLOPE’ ,T3,/, 

* 8X, 'LEVEL* ,17X, 'ST, SL.' ,15X, 'END S L.',/) 

1201 FORMAT (5X,l3,5X,Ai,2X,2(12X,F10.i)) 

1202 FORMAT (16 , ?X, 2( 1 2X ,F10 .D )" " 

1400 FORMATflX, » ASPECT OUTPUT (CLOCKWISE FROM NORTH) » ,914) 
1500 FORMATC PRINT OPT, *',I2,' SLOPE LENGTH PRINT OPT. 
X ' 0 - NO, 1 - YES') - *" ■ — 

2001 FORMAtC 44 ERROR 4* SLOPE *',FI2.I,' IX «',I6,» 

2002 FORMATC *4 ERROR • 44 ELEV. «',F12.1,» IX •',16,' 

2005 FORMATC 1 *44 ELEVATTOir DISTRIBUTION FF4 +-,//, ’ 

* ' LEVEL TOTAL CELLS') 

2006 FORMAT(10X,I3,F15.0) 


',15,' TO ',15,/, 
TO ',15,/) 

' ,F6,t, 


or 


*' , 12 ,/, 


IT 

IT 


'' ,14) 
"»I*> 


C 

c 

c 

c 

c 

c 

c 

c 

c 

-C- 

c 

c 


CARD INPUT - 


CARD 1 
COLS 


1-5 


(15) 


COLS 6-10 


CO L S U-13 
COLS 16-25 
COL 30 


ROW NUMBER ON INPUT FILE TO START 

PROCESSING. DEFAULT « l i 

STARTING COLUMN NUMBER TO PROCESS 
INPUT FILE. 

TWFUT FILET - - 


ON 


(F10.2) CELL SIZE. 

(ID PROCESSING OPTION - 

0 ^ PROCESS SLIT* 

i • ELEVATION DATA ONLT 


7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

-? 

7 

7 

7 

7 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 
46 

46 

47 
‘48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
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60 

61 

62 

63 

64 

65 

66 
67 
60 

69 

70 

71 

72 

73 

74 

75 

76 

77 
70 
79 
00 
01 
02 
03 

84 

85 

86 

07 

08 

89 

90 

91 

92 

93 

94 

95 

96 

97 
90 
99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

it;o 

m 

na 

ni 

114 

115 

lift 

117 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c. 

c 

c. 

c 

c 

c 

c 

c 

c 

c 

0 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CARD 2 
COLS 
COLS 


1-5 

6-10 


COLS 11-15 


2 ■ SLOPE, ASPECT, SL. LENGTH ONLY 


(15) STARTING RON NUMBER FOR OUTPUT. 
(15) ENDING RON NUMBER FOR OUTPUT. 
(15) COLUMN BIAS FOR SHIFTING OUTPUT. 


CARD 3 - (COMPUTES A MAX OF 254 LEVELS) 


COLS 1-8 
COLS 9-16 
COLS 17-24 

CARD 4 
COLS 1-5 


(Ft. 0) 
(F0.O) 
(FR.O) 


(is: 


STARTING ELEVATION LEVEL 
ENDING ELEVATION LEVEL 
ELEV. LEVEL INCREMENT 


NUMBER OF SLOPE LEVELS, INPUT TO FOLLOW. 
MAX 20 


CAROS 5 
COLS 
COLS 
COLS 
COLS 
COLS 
COLS 
COLS 
COLS 
COLS 
COLS 


- (6) 
1-8 
9-16 
17-24 
25-32 
33-40 
41-48 
49-56 
57-64 
65-72 
73-00 


(F8.0) 

(F8.0) 

rrt.oi 

(Ft. 0) 

(F0.O) 

(Ft. 07 
(F8.0) 

(F8.0) 

(Ft .01 

(FB.O) 

INPUT t 
CIO FFR-CARDT 


START SLOPE LEVEL (1) 

START SLOPE LEVEL (2) 

START SLOPE LEVEL (31 
START SLOPE LEVEL (41 
START SLOPE LEVEL 
START' SCOPE” LEVEE 
START SLOPE LEVEL 
START SLOPE LEVEL 
START SCOPE” LEVEL* T9T 
START SCOPE LEVEL (10) 

OR 2 CARDS DEPENDING ON CARD 4. 


(5) 

CS7 

(7) 

(•) 


CARO 7 
COLS 
COLS 
COLS 
COLS 
COLS 
COLS 
COLS 
CQLS 
COLS 


ASPECT OUTPUT PARAMETERS 


1-2 

3-4 

5-6 

7-8 

9-10 

11-12 

IMF 

15-16 

17-18 


crrr 

(12) 

( 12 ) 

TT2T 

(12) 

(12) 

titr 

( 12 ) 

(12) 


— NUMBER COPE F O R” U N D ET E R M INED 
NUMBER CODE FOR NORTH 
NUMBER CODE FOR NC 

NtmSER CODE" rOR-E — 

NUMBER CODE FOR SE 
NUMBER COaE FOR 8 

NUMBER COOE FOR N 
NUMBER CODE FOR MW 


CARD 8 - PRINT OPTIONS 

0 a NO 1 « YES 

COL 2 PRINT ELEVATION, SLOPE AMDA1PECT NAPS 
COL 4 PRINT SLOPE LENGTH CALCULATIONS. 


CALL CSINIT 
CALL PAGES 

r.imrn 

READ (4, 1002) IRON , LROW, ICBIA8 
IF (IO PT.EO.2) GO tO 5 

ir (IOPT.EO.l) GO TO 10 
READ(4 , 1002) NSL 

HEAD! 4 , r&WT*B»CTTrrm7W*fL ) 

READC4, 1004)IA8P 


7 

60 

7 

61 

7 

62 

7 

63 

7 

64 

7 

65 

7 

66 

7 

67 

7 

60 

7 

69 

7 

70 

7 

71 

7 

72 

7 

73 

7 

74 

7 

75 

7 

76 

7 

77 

7 

78 

7 

79 

7 

80 

7 

81 

7 

82 

7 

83 

7 

84 

7 

85 

7 

86 

7 

87 

7 

88 

7 

89 

7 

90 

7 

91 

7 

92 

7 

93 

7 

94 

7 

95 

7 

96 

7 

97 

7 

91 

7 

99 

7 

100 

7 

101 

7 

102 

7 

103 

7 

104 

7 

105. 

7 

106 

7 

107 

7 

ioe 

7 

109 

7 

110 

7 

111 

7 

112 

T 

113 

7 

114 

7 

115 

TUI 

7 

117 
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116 

10 

RF.AD(4,1004) IPO.ISLOP 

7 

118 

119 


IF finPT.F.0.2) GO TO 38 

7 

119 

120 


FNDELV * ENDRLV+KHCir-1. 

7 

120 

*21 

20 

WEI. * C C ENDSLV-STEIiV) /EHCR)42# 

7 

121 

122 


IF fNKL.LE.254) GO TO 30 

7 

122 

123 


FNCR * RMCR+10. 

7 

! 73 

124 


GO TO 20 

7 

124 

125 

30 

XSELft) * 0. 

7 

125 

126 


XEEL(l) * STELV 

7 

126 

127 


XSRL(2) • STELV 

7 

1 . 27 

126 


XFF.f.( 2 ) c XRELC1) ♦ RNCR 

7 

128 

129 


DO 35 Is 3 , NEL 

7 

179 

130 


XSRL(T) s XSRLCI-n+RWCR 

7 

130 

131 


XRFLCtJ = XRRLtl-i ) ♦ FNCR 

7 

131 

132 

35 

CONTINUE *• 

7 

132 

133 


XREL(NFL) ■ 99999. 

7 

133 

134 


IF (inPT.EO.l) GO TO 42 

7 

134 

135 

38 

NL * NSL-1 

7 

135 

>‘36 

. 

DO 40 I * l , HL 

7 

136 

137 


XESL(I) * XSSLCI+1 ) 

7 

137 

13* 

40 

CONTINUE 

7 

136 

130 


XESL(NLTl) » 99999. 

7 

139 

140 

C 


7 

140 

141 

42 

CONTINUE — v ' ■” — 

7 

141 

142 


NROW = LROW-IROW+1 

7 

142 

143 


NCOL * LCOL-ICOL+1 

7 

143 

144 


IF (ISR.EO.03 ISR«1 ' 

7 

144 

145 


NR * NROWtISR-1 

7 

145 

146 


IC « TCOL ♦ ICBI AS 

7 

146 

147 


LC * LCOL ♦ ICBIAS ‘ • ~ * • 

7 

14? 

148 


WRITRC 5.1100) ICOL , LCOL. ISR, CELL* lOPT.NROH# IRON# LROW, NCOL » IC#LC 

7 

148 

149 


IF (IOPT.EG.2) GO TO 150 

7 

149 

150 


WRITE(5 # 1199) ENCR , STELV, XSEL(NEL) " ~ - - 

7 

ISO 

151 


IF (NEL.LE.26) GO TO 145 

7 

151 

152 


DO 140 I s 1 .NEL 

7 

152 

153 


WRITE (5 , 170?) t .XSRLriTrXWLTTT " * ~ 

7 

153 

154 

140 

CONTINUE 

7 

154 

155 

146 

GO TO 150 

7 

155 

1 56 

DO 1^7 X ® 1 f NEL 

. - T 

156 

157 


IL * LSL(Itl92,8) 

7 

157 

158 

147 

WRITE (5.1201) I, IL, XSCL.(J) ,XItL(I) 

7 

158 

I 59 

CONTINUE 

7 

159 

160 

C 


7 

160 

161 

150 

CONTINUE 

7 

161 

162 


IF (IOPT.EO. 1 ) GO TP 1701-— — - 

7 

162 

163 


WRITE ( 5 , 1200) NSL 

7 

163 

164 


NL » NSL 

7 

164 

165 


DO 170 X ^ l f NL 

7 

165 

166 


IL = LSL (1+192,8) 

7 

166 

167 


WRITE (5, 1201) I , IL,X5SL(I) ,XESL(I) 

7 

167 

168 

170 

CONTINUE - - • - 

? 

166 

169 


WRITE ( 5 .1400) IA5P 

7 

169 

i 170 

1705 

WRITE C 5 , 1500) IP0.I8L0F 

7 

170 

171 


JfROW « IROW-r •— — , ™ • 

— ' ■ 7 

171 

172 . 


LO * 1 

7 

172 

173 


ISTP * 0 

7 

173 

174 


KT* K 0 “ " " ■ 

7 174 

; 1 75 


M ■ 0 

7 

175 
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176 



N3 a 0 




7 176 

177 



NCI. = NCOL ♦ 1 




7 177 

178 



DO 171 I « 1,1080 




7 178 

179 



IXFLV(I) a 0 




7 179 

180 



IXSLP(I) a 0 




7 180 

181 



IX4SPCT) * 0 ‘ ’ '• 




7 181 

182 



IXltEN ( I ) * 0 




7 182 

193 

171 

CONTINUE 




7 183 

184 

C 


‘ 




7 184 

1H5 

c 

**« COMPUTE DIST BASED ON CELL SIZE 




7 185 

186 



DO 50 I a 1,8 




7 186 

187 



DISK I) a CELL 




7 187 

198 



IF (MOD (1,2) .EQ. 0) DISTCI) » 1 .4142144CELL 




7 IBS 

189 


50 

continue 




7 189 

190 



IXELVC2) a IC 




7 190 

191 



IXSLPC2) a IC 




7 191 

192 



IXASPC2) a IC 




7 192 

193 



IXI.EN(2) * IC 




7 193 

194 



IXELVC3) = LC 




7 194 

195 



iXSt PC 3) S LC 




7 195 

196 



TXASPC3J ■ LC 




7 196 

197 



IXLENC3) ■ LC 




7 197 

198 

c 






7 198 

) 99 

c 

♦ ** READ INPUT FILE INTO ROTXTWG BUFFER 




7 199 

200 

c 

*** RFPE4T FIRST AND LAST RECORDS 




7 200 

201 

c 






7 201 

202 


100 

CONTINUE 




7 202 

203 



IFCM .EQ. 1) GO TO 110 




7 203 

204 



KTR * KTR+1 




7 204 

205 



IF CKTR.GT.NR) GO TTr TOS ■ * - 




7 205 

206 



CALL VSRDCIFN, 2160, AR, 1ST) 




7 206 

207 



IF (KTR. LT, ISP) GO TO 100 




7 207 

208 



IFCI5T .GT. DJ Gff TITTJO 


- - * - 


7 208 

209 



1ST « -1ST 




7 209 

210 



GO T0C1 10, 110, 105, 103,110,400), 1ST 




7 210 

211 


105 

CONTINUE - * ’ ' 


* 


7 211 

212 



ISTP a 1 




7 212 

213 



LO a 3 




7 213 

214 


110 

CONTINUE " • ---- 

~ *' 


-- 

7 214 

215 



IY * AR(1) ♦ 1 




7 215 

216 



M a M ♦ 1 




7 216 

217 



N3 * N3 M ‘ ~ - 




-3 217 

218 



IFCN3 .GT. 3) N3 * 1 




7 218 

219 



N2 * Ml • 1 




7 219 

220 



XFCN2 .EQ. 0) N2 « 3 ' ‘ 




7 220 

221 



N1 ■ N2 - 1 




7 22! 

222 



irCNl . EQ . 0) N1 a 3 


* 


7 222 

223 



IFCISTP .EQ. II 5UTirT75” ~ ~ 

... 

.... 


7 723 

224 



DO 120 I a l.LCOL 




7 224 

225 



ELEC Cl ,N3) a ARCH-1) 




7 225 

226 


120 

CONTINUE 




7226 

227 



GO TO 145 




7 227 

228 


125 

CONTINUE 




7 231 

229 



DO 13X1 T a r.TTCOL ~ 


'* - • — - 


1 ”229 

230 



ELFCCI.N3) a ELECCI,N2) 




7 230 

231 


130 

CONTINUE 


* 


7 231 

232 


r«S CUWT I HUE 




7 ”212 

233 



IF(M .LI. 3) GO TO 100 




7 2)1 


ORIGINAL PAGE IE 
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234 

235 

236 

237 

238 

239 
740 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 
27t 

279 

280 
281 
282 


C 

*J * 4 

C *** EXTRACT DATA FOR SLOPE C OMP U T A TION 
no 300 I S ICOL,LCOL 
8 * 1-1 

IF (K.LT.ICDL) KSICOL 
L - T ♦ 1 

TF(L , GT. LCOL) L * LCOI. 

XXEL * ELEC(I,N2) 

IF (inPT.EQ.l) GO TO 230 
AP(8) * ELECCK.Nl) 

AP ( 1 ) * ELFC(T,N1) 

APf?) * ELEC(L,N1) 

APC7) s ELEC(K,N2) 

AP(31 * ELECCL.N2) 

AP 1 6) * ELEC(K,N3) 

AP(5) * ELECCI ,N3) 

AP(4D * ELECCL,N3) 

XMAX * 0. 

JJ * 0 

DO 160 J > 1,8 

DIFFCJ) * ASSCXXEL-AP(J)) 

XDIFF * DIFF( J)/0X8T( J) 

IFCXOIFF' .LT. XHAXT~GU~TO LSD 

XMAX * XDIF* 1 
JJ “ J 

• 160 CONTINUE ’ ~ — 

IF (JJ.NE.O) GO TO 180 
SLLEN * 0.0 

GO TD 190 - ’ — - 

180 SLLEN * S0RTCCDIST(JJ)*DIRT(JJ)1 ♦ (HIFFt JJ) 4DIFFI JJ) 5 ) 
190 SLOPE s 30.184XMAX 

IFCXMAX .EO. O.) JJ *' "0 ~~ 

DO 210 K * 1 , N8L 
KK » K 

~IFtSLPPE~.CErXSSti CK) ."ANT I. SL 0PE .L T . TE 8t C*) ) GO TO 220 

210 CONTINUE 

WRITE(5,2001) SLOPE, I, IY 

220 CONTINUE " * — 

IX8LPCN) « KK 
IXASPCN) ■ l A8P( JJ+1 ) 

TXLENtN) = TSET L TN / rO . U0. 5 

IF (IXLEN(N).GT.255) IXLEM(N)*2S5 
TF (inPT.F.0.2) GO TO 290 


C DO ELEVATIONS 

230 CONTINUE 

00 7j ff -n = l ,9Et 

MM * K 


283 IFCXX EL .GE.XS ELCK) .A ND, XXE L .bT.KHL(tt) ) GO TO 270 

285 WRITE(5,2002) XXEL,I,XY 

286 270 CONTINUE 

287 IXELVCN) — 

288 ELF(MM) * ELF(MM)+1 • 

289 290 N s N+l 

2»ff '300 CONTINUE 1 — 

291 NROM « MR0H41 


7 234 
7 235 
7 236 
7 237 
7 238 
7 239 
7 240 
7 241 
7 242 
7 243 
7 244 
7 245 
7 246 
7 247 
7 248 
7 249 
7 250 
7 251 
7 257 
7 253 
7 254 
7 255 
7 256 
7 257 
7 258 
7 259 
7 260 
7 261 
7 262 
7 263 
7 264 
7 265 
7 266 
7 267 
7 268 
7 269 
7 270 
7 271 
7 272 
7 273 
7 274 
7 275 
7 276 
7 277 
7 278 
7 279 
7 .210 
7 281 
7 202 
7 283 
T 284 
7 285 
7 286 
T 287 
7 288 
7 289 
T 290 
7 291 
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292 


TF (MRnW.GT.LROW) GO TO 400 

7 

292 

293 


IXFLVC1) * 4Rf)W 

7 

293 

294 


JXSLP(l) e "ROW 

7 

294 

295 


iXASPfl) * MROW 

7 

,V95 

296 


IXLEN(l) * MROW 

7 

296 

297 


TF C JOPT.EO.2) GO TO 320" * 

7 

297 

298 


CALL VSWR(FFN, 1080, TXELV, 1ST) 

7 

298 

299 


IF ( 1ST , EO. -4) GO TO 350 

7 

299 

300 


IF (IOPT.EQ, 1) GO TO 330 

7 

300 

301 

320 

CALL VSWRfSFN, 1080,1 XSLP# 1ST) 

7 

301 

302 


CALL VCWRUrN, 1080, TXASP, Iff) 

7 

302 

303 


CAfiL V8WR(LFN.10 bO, IXLEN , 1ST) 

7 

303 

304 

330 

M * 2 

7 

304 

305 


Lfl : 2 

7 

305 

306 


IFdSTP .FQ. 0) GO TO 100 

7 

306 

307 


GO TO 400 

7 

307 

308 

350 

CONTINUE 

7 

306 

309 


IF (IOPT.F.O.2) GO TO 360 

9 

309 

310 


CALL VSSR (EFN ,1,1) 

7 

310 

311 


IF (IOPT.EQ. 1) GO TO 400 

7 

311 

312 

360 

CALL V$SR(SFN,l,l) 

7 

312 

313 


CALL V6SRCAFN, 1.1) 

7 

313 

314 


CALL V«SR(LFN,1,1) 

7 

314 

315 

400 CONTINUE " ' " ~ * 

7 

315 

316 


IF (IOPT.EQ. 2) GO TO 405 

7 

316 

317 


CALL VSEF(EFN) 

7 

317 

318 


IF (IOPT.EQ. 1) GO TO 410 

7 

318 

319 

405 

CALL V8EFCSFN) 

7 

319 

320 


CALL V6EF(AFN) 

7 

320 

321 


CALL VSEFTLFN) ‘ ^ " ' 

7 

321 

322 

C 


7 

322 

323 

C 


7 

323 

324 

c 

■ ,*■ " “ *" ** 

7 

324 

325 

. 410 

CONTINUE 

7 

325 

326 


IF (IPO.EQ.O) GO TO 800 

7 

326 

327 


IF CinPT.Eff.TJ GO TCT'590 " * 

7 

327 

328 

c 


7 

328 

329 

c 

*** PRINT ELEVATION DATA *** 

7 

329 

330 


IF CNEL.CE.267 GO TO *W 

7 

330 

331 


NRX * (HCOL<’29)/30 

7 

331 

332 


DO 480 I ■ l , NRX 

7 

332 

333 


CALL V5FTTEFfT|T(r,TT ~~ — ‘ 

7 

333 

334 


CALL PAGE* 

7 

334 

335 


IS * (I-l)*30 ♦ ICOL 

7 

335 

336 

■* “* • 

IE ■ IS+29 - - 

7 

336 

337 


IP CIE.GT.LCOL) IE*LCOL 

7 

337 

338 


PARM ( 1 ) ■ ELEC 1 ) 

7 

338 

339 


"Piwrro TS-ELEC2) 

7 

33.9 

340 


ISA = IS+IC8I AS 

7 

340 

341 


IEA * IE+IC8IAS 

7 

341 

342 


ViRITE(5,5I0)PARW(n,PAWRC2),II»;iEI - - 

7 

342 

343 


IS • IS-IC0L44 

7 

343 

344 


IE ■ IE-ICOL+4 

7 

344 


«5’ 

346 

347 

«r 

349 


JF(IST.EQ.-3 .OR. IST.E0.-4) GO TO *90 
IY ■ IXELVCl) 

WRITE! svswrn ; i 
GO TO 420 


345 

346 

347 

341 

349 
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350 480 CONTTMUP 7 350 

351 MRX = (NCDI.+ 1 19J/120 7 35 t 

352 GO TO 590 . 7 352 

3^3 C 7 753 

354 490 NR X s (NCOL+1 19J/120 7 354 

355 DO 500 I = 1 , NRX * 7 355 

356 CALL VSFCCEFN.O#!) 7 356 

357 CALL PAGE* 7 357 

358 TS * fT-l)M20tICOL 7 358 

359 IE * IS ♦ 119 7 759 

360 TFCIE .GT. LCDL) IE « LCOL 7 360 

361 PAPH(I) = ELEC 1 ) 7 361 

362 PAPM(2) * FLRC2) 7 362 

363 ISA = I5+XC4IAS 7 363 

364 TEA * IE+ICBIAS 7 364 

365 WRITE(5,510) PARM(1),PARM(2) , ISA, IE A 7 365 

366 IS * IS-ICOL ♦ 4 7 366 

367 IE * IE«ICOL +4 7 367 

368 520 CALL VSRDtEFN, 1 080 , IXELV , 1ST) 7 368 

369 IFCIST.E0.-3 .OR. TST.E0.-4) GO TO 500 7 369 

370 IY « IXELVC1) 7 370 

371 DO 525 F * IS, IE 7 371 

372 525 IXELVCK) « LSLCIXELV(K>4l92,8) 7 372 

373 WRITE(5,530) IY, (IXELVTK) ,K*TS,TE7 7 373 

374 GO TO 520 7 374 

375 500 CONTINUE 7 375 

376 C .. . 7 376 

377 C •** PRINT SLOPE DATA *MM******M*M»»»**9****9*******M*f **»* 7 377 

378 C 7 378 

379 590 CONTINUE - - 7 37*3 

380 IF (lOPT.EG.l) GO TO 999 7 380 

381 NRX * (NCOL+119J/120 7 381 

382 DO 600 I ■ l, NRX * 7 382 

383 CALL VSFC(SFN,0,1) 7 383 

384 CALL PAGE* 7 384 

385 IS = CI«im7U + TCmr“*' - 7 38S 

386 IE = IS + 119 7 386 

387 IFCTE .GT. LCOL) IE s LCOL 7 3R7 

389 PARK ( l ) * SLOtlT 7 388 

389 PARMC2-) * SLOC2) 7 389 

390 ISA c IS+ICBIAS 7 390 

391 IEA a IE+TCBTAS” — “ * ~ 7 391 

392 WRITE! 5 ,510) PARK ( 1 ) , PARM (2) , ISA, IEA 7 392 

393 IS * IS - ICOL ♦ 4 7 393 

394 IE *- IE - ICOIT ^ 6 — ~ 7 394 

395 620 CALL V*RD (SFN , 1080, IXSLP, 1ST) 7 395 

396 irCIST.EG.-3 ,OR. XST.EO.-4) GO TO 600 7 396 

397 IT S rXSLPtn • 7 397 

398 • DO 625 K * IS, IE 7 398 

399 625 IXSLPCK) * LSLCIXSLPC K)+l92,8) 7 399 

400 NITITrt5,S307 IT, ClXSLPTI),K»IS,IE) 7 400 

401 GO TO 620 7 401 

< 402 600 CONTINUE ' 7 402 

403 — C — — — — 7 403 

404 C *** PRINT ASPECT ANGLE DATA ************************************* 7 404 

405 C 7 405 

• 40« TTO~700 X * lYWRX " 'T 498 

, 407 CALL V*FCCAFN,0,1) 7 407 
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P*GE 17 01/04/78 47133000 VORTXTI FTN IV(G) D148A1 0927 HOURS 


400 


CALL PAGES 


7 

4flP 

409 


IS * (1-11*170 + I COL 


7 

409 

410 


TE * IS 4 U9 .. . . 


7 

41 r 

411 


IFCIE .GT. LCOL) IE « LCOL 


7 

41 1 

412 


PARM(l) * ASP ( 1 ) 


7 

4 1 i 

413 


P ARM ( 2 ) a ASP(2) 


7 

4 1 i 

414 


ISA a IS+IC4rAS 


7 

4 l 4 

415 


IEA = IE+IC8I AS 


7 

415 

416 


WRITE(5 # 51 0 ) PARM(l) ,PARM(2) , ISA, IEA 


7 

4 1 «• 

417 


IS a IS - ICOL ♦ 4 


7 

417 

418 


IE = IE - ICOL ♦ 4 


7 

418 

419 

720 

CALL VSRDCAPN, t080,IXASP,fSTJ 


7 

419 

420 


lF(IST.E0.-3 .OR. 1ST. EG. -4) GO TO 700 


7 

420 

421 


TY a IXASPCl) 


7 

421 

422 


WRITEC5, 531) TY,(IXASP(K),K*IS,IE) 


7 

422 

423 


GO TO 720 


7 

42; 

424 

700 

CONTINUE 


7 

424 

425 

800 

CONTINUE " ' 


7 

425 

426 


IF (ISLOP.E1.0) GO TO 999 


7 

426 

427 

C 

♦ 


7 

4 2' 

428 

C *44 PRINT SLOPE LENGTH »44664**4*4*46«46666 


7 

428 

429 

C 

30 COLS PER PAGE 


7 

429 

4 30 

C 



? 

4 3 ^ 

431 


NRX a f NCDL+29) /JO 


7 

4 1 i 

432 


DO 900 I a 1 r NRX 


7 

4 ’ 2 

433 


CALL V6FC ( LFN ,0,1) 


7 

4 Ji 

434 

- 

CALL PAGE! ' 


7 

4 3 •» 

435 


IS « (1-1)430 + ICOL 


7 

435 

436 


IE « IS429 


7 

436 

437 


IF CTF.GT.LCOL) IT * LCOIT “ ~ 


7 

4 3* 5 

438 


PARM(l) = SLEN(l) 


7 

43* 

439 


PARMC2) a SLENC2) 


7 

439 

440 


ISA * IS+TC9IA S * 


7 

4 £0 

441 


IEA * micaiAs 


7 

441 

442 


WRITE(5,510) PARMCl) ,PARM(2) ,ISA,IEA 


7 

442 

443 


WRTTET5 , 511 T 


7 

443 

444 


IS * IS-ICOL+4 


7 

444 

445 


IE a TE-ICOL+4 


7 

445 

446 

920 

CALL VtRBCLFN, 10907TXTjEN,IBT) 7 ” 


7 

446 

447 


IF (ISJ.E0.-3 .ON. I8T.E0.-4) GO TO 900 


7 

447 

448 


IY * IXLENC1) 


7 

446 

449 


WRITE (5 ,540 ) IYrTTXL^TffCIO » K*T5 , IE1 


7 

449 

450 


GO TO 920 


7 

4 50 

451 

900 

CONTINUE 


7 

451 

452 

999 

CONTINUE 


7 

•452 

453 

C 



7 

453 

454 

C 

PRINT ELEVATION PREOUENCT DIITRH. 


7 

454 

45T 

~cr ■ 



7 

45* ! 

456 


IF (I0PT.E0.2) GO TO 1000 


7 

456 

457 


WRITE (5.2005) 


7 

457 

Iff 

mr 9io i • i « nec' 


7 

450 

459 


MRXTE(S,2006) X,ELP(I1 


7 

459 

460 

914 

CONTINUE 


7 

460 

461 

1000 

STOP 


7 

461 

462 


END 


7 

462 


ENTRY/COMMON BLOCK NAMES 

_ r oooooo TDPoini 

EXTERNAL NANE8 
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APPENDIX D 


GENERAL INFORMATION CONCERNING DATA SOURCE 


GENERAL INFORMATION CONCERNING DATA SOURCE 


The NCIC, Office of Research and Standards, Reston, Virginia 
and the DMATC, Washington, D.C. were contacted in regard to the 
digital terrain napes . It was learned that the DMATC prepares 
two 1° x 1° matrices for each 1 : 250 ,000-scale quadrangle map. 

The x, y sheet corners and all subsequent readings of each half 
of a quad map are referenced to an arbitrary coordinate system 
(ACS). ‘The origin of this ACS is an arbitrary reference point 
that is different for each 1° x 1° area. The approach used at 
the ERL is to reference the (x,y) plate readings directly to 
the UTM grid and then compute coordinates of points in the UTM 
system directly from the readings. Thus, an angle correction 
is required for each 1° x 1° area. 

Further discussions revealed that a photographic process 

is used to produce a map x from the 1 : 250 , 000-scale series whereby 

a "contour" line has a "ditch" type depression. The operator of 

the digital graphic recorder follows this ditch in the actual 

extraction of contour data. Further processing involves inter- 

% 

polation of a planar nature. Error sources in the DMATC process 
are resolution (<!0 mils) and operator accuracy (less than 
resolution error) , 

It was also verified from the DMATC that there is an "edge" 
matching problem on the right and left of all quad maps and 
that the data exhibits a "saw-tooth" effect at the top or bottom 
edges of quad maps. Since the actual extraction of planar data. 



(x-y plate coordinate readings in inches) is from a non-square 
map (latitude-longitude map) , then the subsequent transformation 
of these planar data to the regular shaped UTM coordinate system 
results in these edge and border problems. These problems have 
been encountered in the research efforts at the ERL. It should 
be pointed out that these problems are not encountered if the 
subsequent analysis programs (such as data base programs) are 
compatible with the latitude- longitude coordinate system. 



